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ABSTRACT 



i-£h ■ It has been shown that the dissipative gas infall during galaxy formation has 

I*] the capability to modify the shape of dark matter halos. In this paper we perform 

Q . the first detailed analysis of particle orbits in a cosmological dark matter halos to 

i—{ • understand how and why baryons alter its shape. We perform a series of numerical 

c/3 \ experiments where we grow a baryonic core inside a live dark matter halo extracted 

from a cosmological simulation. We follow the orbits of more than 300 particles with 
more that 50000 timesteps. Our results clearly show that the dissipational component 
, is responsible for repeatedly deflecting orbits which visit often the center of the system. 

J> ■ Moreover the gravitational potential time dependence associated with the growth of 

OO ' the baryonic mass, shifts the frequencies of the orbits, making them extremely chaotic. 

t N *»- . This randomization makes the orbits explore a large phase space. When this effect 

■ takes place for a significant number of orbits it will be manifested in the density 

distribution as an approach to a rounder configuration. As a consequence, the influence 

t^J- i of the central mass on the shape of the phase space decreases with increasing distance 

f^S ' from the center. We discuss the importance of future analysis of controlled experiments 

t***'- \ (i.e. using analytic potentials instead of live DM halos) to better decipher the dynamics 

■ of this phenomenon. 
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1 INTRODUCTION 

Currently, the favored cold dark matter (CDM) model of 
structure formation suggests that dark matter (DM) halos 
of galaxies and clusters should admit triaxial density pro- 
file (Barnes & Efstathiou 1987, Dubinski & Carlberg 1991, 
Warren et al. 1992, Jing & Suto 2002, Moore et al. 2004). 
This expectation is well-accepted, but (in principle) it could 
also be tested against observations. In this respect, triax- 
ial dark halos surrounding both elliptical and disk galaxies 
could, not only provide clues about the underlying physics 
of galaxy formation, but also impose constraints on differ- 
ent formation scenarios (Dubinsky & Carlberg 1991). On 
the other hand, the recently observed coherence of the tidal 
stream of the Sagittarius dwarf spheroidal suggests a nearly 
spherical halo configuration for the Milky Way (Ibata et al. 
2001, Mayer et al. 2002, Majewski et al. 2004), in contradic- 
tion with what has been obtained by cosmological simula- 



tions (i.e. Jing & Suto 2002, Allgood et. al 2005, Maccio et 
al 2007). 

This inconsistency could possibly be resolved by care- 
fully examining the effect gas cooling has on halo shapes. In 
1994 Dubinsky studied the effects of dissipative gas infall on 
the shape of DM halos by adiabatically growing central-mass 
concentrations inside initially triaxial DM halos. A more ex- 
tensive study was made recently by Kazantzidis et al. (2004; 
K04 hereafter), where a number of cosmological gas dynam- 
ics simulations were performed; the main result was that ha- 
los formed in simulations with gas cooling are significantly 
more spherical than corresponding halos formed in adiabatic 
simulations. 

Aiming to extend the aforementioned investigations 
(Dubinsky 1994 and K04), this paper attempts to address 
the question of why and how the gas dissipation alters the 
halos shape. The main addition is that we approach the 
problem from a dynamical point of view. 
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The relation between particle orbital properties and the 
shape of bound structures has been already studied by sev- 
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oral authors (i.e. Jesseit, Naab & Burkert 2005, Naab, Jesseit 
& Burkert 2006, Valluri et al 2007). In this paper we present 
the first broad analysis of particle orbits in cosmological DM 
halos, and we investigate how these orbits modify their be- 
haviours when a baryonic core is grown inside the original 
triaxial dark matter halo. 

The analysis of the orbits is two-leveled. Firstly, we 
study them in a statistical fashion, investigating quantities 
like the evolution of axes ratios of the configuration. Sec- 
ondly, we examine the behaviour of individual orbits; we ask 
how their shapes are affected when there is a dissipational 
component at the center of the configuration, but we also 
employ a frequency analysis scheme to quantify the chaotic- 
ities, and we attempt to investigate how chaos connects to 
the change of their shapes. 

This paper is organized as follows: in section 2 we 
present our codes and modeling; results are presented in 
section 3; section 4 contains our conclusions and description 
of follow-up for this study. 



2 NUMERICAL SIMULATIONS 

All simulations were performed using PKDGRAV, a tree- 
code for cosmological simulations of N bodies, written by J. 
Stadel and T. Quinn (Stadel 2001). This code uses spline 
kernel softening, for which the forces become completely 
newtonian at two softening lengths. Individual time-steps 
for each particle are chosen to be proportional to the square 
root of t he so ftening length e, over the acceleration, a: 
Ati — rjy/e/al. Forces are computed using terms up to hex- 
adecapole order and a node-opening angle 8 which we set to 
0.6. 

The original dark matter halo was extracted from an 
existing iV-body simulation where we had followed the evo- 
lution of 250 3 particles in a volume of 20 Mpc with a soft- 
ening e = 0.6kpc. We adopted a flat ACDM cosmology with 
parameters from the first year WMAP results (Spergel et 
al. 2003): matter density Qm = 0.268, baryon density fib = 
0.044, Hubble constant h = # /(100km s _1 Mpc _1 ) = 0.71, 
and a scale-invariant, Harrison-Zeldovich power-spectrum 
with normalization erg = 0.9. A detailed analysis of this sim- 
ulation has been already presented in Maccio et al. (2007). 

At redshift zero we selected a galaxy-like halo having 
the following characteristics: M vir = 1.7 x 1O 12 M0, R vir = 
248 kpc and about 110000 particles within the virial radius. 
This halo is well-fitted by an NFW (Navarro, Frenk and 
White 1997) profile with a concentration parameter of 9.5. 
We also computed the shape of this halo with the same 
technique described in Maccio et al. (2007), getting for the 
axis ratio b/a = 0.82 and c/a = 0.71, where a,b,c are the 
long intermediate and short axis respectively. For our further 
studies we cut from the original simulation a cubic region 
of lxlxl Mpc h _1 centered around the selected halo. This 
region sets the new simulation volume for the baryonic mass 
growth experiment. 

2.1 Baryonic core modeling 

The investigation of galaxy formation inside dark matter ha- 
los is a challenging problem. In hierarchical models a collec- 
tion of smaller scale fluctuations (protogalaxies) merge into 



forming a smooth, centrally concentrated galaxy. Both grav- 
itational collapse and dissipation occur at the same time; if 
the cooling time of dissipation is shorter than the gravita- 
tional free-fall time, then gas will sink into the center of 
the system to form a luminous galaxy. Following Dubinski 
(1994) we modeled the formation of the baryonic core of our 
galaxy-model by slowly growing a mass distribution at the 
center of a previously selected dark matter halo. The mass- 
growth follows a linear law, until it reaches a maximum mass 
comparable to the expected luminous mass of galaxies. We 
chose Mb — 1.8 x 10 10 M Q for the final baryonic mass and 
T = 7 Gyrs as growth time (this is the time elapsed from 
z — 1 to z — in our ACDM universe). We also performed 
two more simulations with a smaller final mass for the bary- 
onic core (keeping constant the growth time), namely we 
used a 1/4 and 1/16 of the previous value for Mb- For com- 
parison purposes, we also performed a simulation with no 
growing mass at the center, where we simply allowed our 
halo to evolve in isolation (i.e. once extracted from the N- 
body simulation) for 7 Gyrs. After that time we stopped the 
mass growth on the baryonic core, and we let all the afore- 
mentioned simulations evolve for an additional 7 Gyrs, in 
order for our halos to reach a well-relaxed configuration. 

In first approximation, we treated the luminous compo- 
nent of the galaxy as a single particle within the simulation. 
At time T = we gave to this new particle the same coor- 
dinates and velocities components of the dark halo center of 
mass. This particle was then free to move in the usual self- 
consistent fashion of an iV-body simulation, but, in practice, 
it quickly settled at the center of mass of the dark matter 
halo because of dynamical friction. 

One may object that this zeroth-order approximation 
might generate unphysical effects; for example having only 
one particle at the center of the configuration may cause 
unphysical close encounters with the dark matter particles, 
with subsequent scattering of the last. To assure that this 
was not the case in our experiments, we quantified the im- 
portance of these effects in order to separate physical phe- 
nomena from numerical artifacts generated by the approxi- 
mations of our schemes. 

To this end we performed two additional simulations. 
Keeping constant the growth-rate of the central mass we: 
i) increased the kernel softening for the growing particle to 
four times the one on the surrounding dm particles, from 
0.6 kpc to 2.4 kpc, and ii) instead of using one particle 
we split its mass into 500 particles. The initial spacial size 
of this particle-cluster was 5 kpc and all the particles in- 
creased their masses following the same growth-rate. As is 
clearly stated in the discussion section our results suggest 
that there is no significant, quantitative or qualitative, vari- 
ation between these two simulations and the original one 
(one particle and small softening) . This encourages us to be 
confident about the physical (as opposed to the numerical) 
origin of our results. 



3 NUMERICAL EXPERIMENTS 

We performed a number of simulations of collapsing halos 
with and without a growing central baryonic mass. The den- 
sity configuration of the initial conditions was triaxial with 
axes a > b > c. In these simulations we selected the same 
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sample of 300 particles to analyze. For these particles we 
recorded their positions and velocities for more than 55.000 
time-steps throughout the total evolution time (14 Gyrs). 
The recording time-step of the orbits was set equal to the 
smallest time-step achieved by each particle during the force 
integration. Those tracked particles uniformly cover the en- 
ergy distribution of all the particles at T = and their 
distance from the center ranges from 4 to 150 kpc. We also 
tracked the orbits of all the growing particles although we 
did not use them for the orbits modification analysis. 

We investigated carefully both the behaviours of indi- 
vidual orbits, and the way these behaviours affect the mor- 
phological nature of the central regions of our systems. The 
evolutions of the axes ratios of both dissipational and dissi- 
pationless regimes were determined by employing an analy- 
sis of their inertial tensors. 

The general phenomenology seen already in previous 
works (K04) , emerged in our simulations too; in the absence 
of a central baryonic mass, the original triaxiality insists 
throughout the evolution of the system. On the other hand, 
when the central dissipational component is present, there 
is a rapid increase of the values of the axes ratios b/a and 
c/a which is more evident at the central parts of the sys- 
tem, and intensifies as time evolves (Fig. 1). This increase 
is manifested as a visual "rounding" of the contours of the 
density configuration (Fig. 2). 

Several investigators (Norman et al., 1985, Gerhard and 
Binney, 1985, Barnes and Hernquist, 1996,Merritt and Quin- 
lan, 1998, Valluri, and Merritt, 1998) have provided a dy- 
namical explanation for a similar behaviour occurring in the 
context of galactic systems with central black holes. They 
have shown that black holes alter the potential at the center 
of the system to a rounder one, and they modify several of 
the box orbits and tube orbits with low angular momentum 
(these orbits naturally visit the center of the system and 
serve as important building blocks for the triaxial backbone 
of ellipticals.) Also, the phase space at the center should 
follow the potential and become rounder too. Since a signif- 
icant number of orbits behave chaotically they will attempt 
to cover all the phase space energetically available to them. 
Therefore one expects the motion of orbits to build density 
configurations with rounder shapes. 

A comment about timescales of chaotic diffusion may 
be useful here. In general, chaotic orbits attempt to access 
all the phase space energetically available to them. However, 
the timescale necessary to achieve that, relates closely to the 
dynamical structure of the phase space. When parts of the 
phase space are occupied by regular orbits, chaotic orbits 
located very close to them try to behave like regular too, 
and remain constrained in restricted regions of the phase 
space for a long time, often longer than the age of the sys- 
tem (Contopoulos 2002, Efthymiopoulos et al. 1997, Sideris 
2006); these orbits are usually called "sticky" or "weakly 
chaotic". Sticky orbits may need a significantly long time 
to escape into a wider chaotic sea and thus become "wildly 
chaotic" . Only then they can occupy a much broader phase 
space area. On top of that, in a three-dimensional phase 
space although different chaotic areas are interconnected, it 
is the location and extent of the regular regions which de- 
termines how fast the diffusion of chaotic orbits occurs, and 
whether different parts of a chaotic sea can work as "bottle- 
necks" of the evolution. 



Another issue is the influence of the time-dependence. 
In the experiments we performed there was the significant 
time-dependence associated with the growth of the central 
mass. In time-dependent regimes every energy at every time 
is associated with its own phase space structure. Orbits do 
not conserve energy, they move from energy level to energy 
level, thus they visit a plethora of different phase spaces; this 
means that their nature can actually alternate from regular 
to chaotic and viceversa (Kandrup et al. 2003). 

We investigated the chaoticity of a sample of 300 or- 
bits randomly chosen from the orbits located at the central 
part of our density configuration. The maximum initial ra- 
dius of the sample was about 150 kpc. The first question 
we asked was: how does the chaoticity of the dissipational 
system compare to the chaoticity of the dissipationless one? 

Regular orbits are associated with a limited number of 
frequencies of motion (these being the main frequency per 
degree of freedom and its harmonics). Theoretically these 
frequencies should be singular, e.g. delta-functions. On the 
other hand, chaotic orbits are associated with many frequen- 
cies; their spectra are broad and theoretically continuous 
(Lichtenberg & Lieberman 1992). 

Numerically one does not have the luxury to integrate 
orbits for infinite or even sufficiently long time. This causes 
some unavoidable numerical artifacts; frequencies of regu- 
lar orbits, although very localized, appear to have small 
but measurable tails in the power spectrum. Frequencies 
of chaotic orbits are numerous but discrete. Despite these 
difficulties, one can usually distinguish easily between reg- 
ular and chaotic orbits, provided that the orbit has been 
integrated for an adequate number of orbital periods. As 
the number of available orbital periods decreases, our con- 
fidence regarding the characterization of the orbit decreases 
accordingly. A method of computing localized frequencies 
in localized time would be invaluable, especially for limited- 
time evolutions, but such method has not been clearly for- 
mulated yet. 

In this spirit, we quantified chaos using a straightfor- 
ward numerical technique associated with Fourier analysis ( 
Press et al. 1993; Kandrup et al. 1997). For every orbit we 
computed its three Fourier spectra, one per degree of free- 
dom. Then (after we sorted the frequencies in descending 
order of their power) we simply counted how many of the 
emerging frequencies consist the 90% of the power spectrum 
(starting with the strongest frequency and moving down to 
weaker ones). (Similar experiments with 80% and 95% gave 
similar results qualitatively) . Finally we added the numbers 
of all three dimensions: 

no. 90 = n x oso + n y o.9o + n z o.90 (1) 

The results of this analysis suggest that, when there is 
dissipation, chaos increases significantly at least close to the 
center of the system (Fig. 3) . The characterization results for 
orbits moving close to the center are reliable since the num- 
ber of orbital periods is about 30-50. As one moves away 
from the center, the number of orbital periods decreases 
and the reliability of the characterization unavoidably de- 
teriorates. Still the trend close to the center is obvious; the 
dissipational component makes the system more chaotic. 

Naturally, one may ask; what does really cause this in- 
crease of orbital chaoticity? And is this increase associated 
closely to the rounding of the shape of the halo at the cen- 
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tral parts of the configuration? To answer these questions a 
number of elements have to be examined. Firstly, one can 
look at the Poincare section of the orbits. Such a section, 
recorded when the orbits cross the v x = plane, reveals 
that when the dissipational component is present, velocities 
increase close to the center (Fig. 4). This is not a surprise 
and is clearly associated with close encounters of the orbits 
with the central baryonic mass; the particles approach close 
to the baryonic mass and get deflected. 

It is reasonable to assume that the close encounters with 
the central mass, and the time-dependence caused by the 
growth of the central mass can combine to cause a shift 
to the frequencies of motion of individual orbits with subse- 
quent randomization of the direction of the orbital evolution. 
This way any regularity or sense of shape disappear, the 
orbit attempts to occupy phase space in many different di- 
rections, thus becoming highly chaotic. This process should 
be manifested as a broadening of the power spectra of or- 
bits (since new frequencies appear fast). A look at the three 
frequency spectra of same orbit moving in the two different 
configurations (with or without a growing central baryonic 
element) makes obvious that the spectrum is much broader 
in the second case (Fig. 5). The existence of the broader 
spectra obviously correlates to the increase of chaoticity ob- 
served in Fig 3. 

Although we have important hints about the dynami- 
cal process followed, one needs a careful quantification of if 
and how the direction of the emerging orbital components 
is really randomized. This effect can be visualized clearly 
for an orbit passing very close to the center. One may ask 
how the distribution of the angles 8 — tan' 1 (x/y) changes 
between the dissipational and dissipationless regime. In a tri- 
axial dissipationless regime such a distribution should have 
some preferences, because there are no intense phenomena 
like strong deflections. On the other hand a randomization 
of the orbit should be associated or approach a distribution 
closer to normal, since all possible angles must have equal 
probability to be visited. Indeed by overplotting the 6 dis- 
tributions for the same orbit in both regimes it is obvious 
that what was expected actually happens Fig 6. 

This randomization will not happen for every orbit in 
the system, but a significant percent of orbits in our sample 
followed the aforementioned process. By visual inspection 
of a number of orbits one can easily identify a significant 
percentage which try to mimic and are restrained in the 
box-orbits fashion. The close encounters with the central 
dissipational components will make these orbits rounder. It 
has to be emphasized that although box-like orbits with big- 
ger amplitudes can in principle experience similar dynamics, 
one has to remember that as the amplitude increases the 
randomization effect decreases simply because the number 
of orbital periods is smaller. This means that the opportuni- 
ties for close encounters with the central gas are decreased. 
This seems to explain why the effect of rounding is limited 
only in the central parts. A second reason should be that 
the effect the central mass has on the rounding of the phase 
space weakens in bigger distances from the center. Box-like 
weakly chaotic orbits with big amplitude do not have enough 
time to "break" and explore many new directions during 
their evolution. In an imaginary system which evolves for 
far longer times we would probably observe a rounding of 
the density in bigger radii too. It is the evolution time of 



the system that puts a limit in the maximum radius it gets 
rounder because of deflections by a central growing baryonic 
mass. 

Figure 1 also shows that the shape modification is not 
constant with respect to time. It is faster in the first 5 Gyrs 
for both the axis ratios: b/a and c/a. This is related by the 
total fraction of box-like orbits within the central region of 
the halo. The 3-dimensional distribution of the mass is inti- 
mately connected to its orbital content. Jesseit et al. (2005) 
showed correlations between orbital fractions and the shape 
parameters as determined by the axis ratios of the inertial 
tensor. To give a rough estimate on the box orbit fraction 
in our initial halo, we selected the merger remnant from the 
original sample of Naab & Burkert (2003) which resembles 
the analysed DM halo the most. Hence we conclude that 
the DM Halo has a box orbit fraction of at least 0.67 for the 
inner 10% , 0.49 for the inner 20 %, 0.39 for the inner 30% 
, 0.33 for the inner 40% , 0.27 for the inner 50% and 0.23 
at 60% of all particles. The fast change in the shape in the 
first Gyrs is related to the high number of box-like orbits 
(i.e. highly modificable orbits) within the central region of 
the halo. 

Some additional experiments were necessary in order to 
corroborate the understanding of the physics. First we per- 
formed the same simulation using a much bigger softening 
parameter e = 2.4. This simulation was necessary in order to 
show that the effects we observed where not caused by the 
numerics of the experiment. Comparison of the evolution of 
the axes ratios show clearly that the larger epsilon does not 
alter the results (Fig. 7). Moreover, the chaoticities of the 
orbits in this regime do not differ significantly from the ones 
evolving in a central growing mass regime (Fig. 9). This is a 
strong indication that the observed effects are not numerical 
artifacts. 

The second experiment was a more realistic version of 
the first one. Specifically, we replaced the central baryonic 
mass with a set of 500 smaller baryonic masses. The to- 
tal mass and the rate of mass growth were kept the same. 
The 500 baryonic masses were localized at the center. By 
comparing the effects of the two dissipational simulations 
one can establish that the crucial aspects of dynamics are 
similar in both cases. Indeed, the evolution of the axes ra- 
tios are similar although they start differing as we approach 
very close to the center (Fig 8) . This is reasonable since the 
deflections are not very strong anymore, therefore the effect 
needs more time to become evident. An analysis of the chaos 
of the system suggests that the chaoticity of the 300 orbits 
we analyzed are in relatively lower level than when we had 
one central mass but still higher than the dissipationless case 
(Fig 9). 

It is very instructive to see how the same orbit changes 
in the four different cases, (a) dissipationless system, (b) 
dissipational system, one central mass, (c) dissipational sys- 
tem one central mass e — 2.4, (d)dissipational system 500 
central masses. Fig. 10 provides the three projections x — y, 
y — z, and x — z, of the same orbit for all the above cases. 
We chose this specific orbit because the difference between 
the dissipationless and dissipational cases are dramatic and 
thus make the phenomenon clearer. The effect we described 
is particularly obvious in the y — z projection. In the dissi- 
pationless case the orbit has a strong radial component and 
passes very close to the center. When there is dissipation 
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the central mass makes the potential rounder in the central 
areas where its influence is strong and it also deflects the 
orbit to several different directions. When e = 2.4 the orbit 
also becomes thicker. For the case with 500 at the center the 
rounding effect is also present. 



4 DISCUSSION AND CONCLUSIONS 

This investigation suggests that the dissipational component 
is responsible for repeatedly deflecting orbits which visit of- 
ten the center of the system. Also the time-dependence as- 
sociated with the growth of the central mass, (probably in 
combination with general nonlinearities of the triaxial con- 
figuration evolution), seem to shift the frequencies of the 
orbits, making them extremely chaotic. These effects ran- 
domize the directions achieved during the evolution of the 
orbit. This randomization makes orbits explore many more 
angles, and therefore they occupy rounder shapes. When this 
effect takes place for a significant number of orbits it will be 
manifested in the density distribution as an approach to a 
rounder configuration. 

This approach can easily serve as an explanation of why 
the rounding is limited only at the central parts of the con- 
figuration. Orbits with extensive amplitude cannot visit the 
center many times during the life-time of the system, thus 
they cannot be deflected repeatedly. This way their direc- 
tions cannot be fully randomized in the sense this happens 
for orbits with smaller amplitude and they never get the op- 
portunity to explore a much rounder space. Moreover, the 
influence of the central mass on the shape of the phase space 
decreases with increasing distance from the center. Since the 
shape of the orbits with bigger amplitude is not significantly 
altered by the existence of the central mass, there is no rea- 
son to expect any morphological alteration on the distribu- 
tion in larger radii. 

It has to be noticed that this interpretation does not 
disagree with previous ones. It simply incorporates into the 
picture the role of extra influences into the evolutionary his- 
tory of orbits by repeated central deflections, as well as the 
time-dependence of the central baryonic mass. Two main 
questions are to be explained in future work in more con- 
trolled experiments: (a) What is the influence of each one 
of the two components? What is the exact role of the time- 
dependence of the central mass, and how does it combine 
with (or affected by) the nonlinearities of the self-consistent 
evolution in a triaxial regime? (b) Is there a general pattern 
in the way the frequencies of motion emerge for individual 
orbits? If yes how does this connect with the average loca- 
tion of the orbit as well as the rate of growth of the central 
mass? The answers to these questions can probably decipher 
the fundamental dynamics of this phenomenon and provide 
important clues for the treatment of the problem using an- 
alytics. 



of Zurich. AVM thanks Doug Potter for his fundamental help 
with all the zBox2 issues. 



REFERENCES 

Allgood, B., Flores, R. A., Primack, J. R., Kravtsov, A. V., Wech- 
sler, R. H., Faltenbacher, A., & Bullock, J. S. 2006, MNRAS, 
367, 1781 

Barnes, J. & Efstathiou, G. 1987, ApJ, 319, 575 
Barnes, J. E., & Hernquist, L. 1996, ApJ, 471, 115 
Contopoulos, G. Order and Chaos in Dynamical Astronomy 

(Springer- Verlag, Berlin, 2002). 
Dubinski, J. 1994, ApJ, 431, 617 
Dubinski, J. & Carlberg, R. G. 1991, ApJ, 378, 496 
Gerhard, O. E., & Binney, J. 1985, MNRAS, 216, 467 
Efthymiopoulos, C, Contopoulos, G. Voglis,N., and Dvorak, R. 

1997, Journal of Physics A 30, 8167. 
Jesseit, R., Naab, T., & Burkcrt, A. 2005, MNRAS, 360, 1185 
Jing, Y. P. & Suto, Y. 2002, ApJ, 574, 538 

Kandrup, H. E., Eckstein, B. L., & Bradley, B. O. 1997, A&A, 
320, 65 

Kandrup, H. E., Vass, I. M., & Sideris, I. V. 2003, MNRAS, 341, 
927 

Kazantzidis, S., Kravtsov, A. V., Zentner, A. R., Allgood, B., 
Nagai, D., & Moore, B. 2004, ApJL, 611, L73 

Lichtcnberg, A. J. and Lieberman, M. A., Regular and Chaotic 
Dynamics (Springer- Verlag, New York, 1992). 

Ibata, R., Lewis, G. F., Irwin, M., Tottcn, E., & Quinn, T. 2001, 
ApJ, 551, 294 

Maccio, A. V., Dutton, A. A., van den Bosch, F. C, Moore, 

B. , Potter, D., & Stadel, J. 2007, MNRAS in press, astro- 
ph/0608157 

Majewski, S. R., Skrutskic, M. F., Weinberg, M. D., & Ostheimer, 

J. C. 2003, ApJ, 599, 1082 
Mayer, L., Moore, B., Quinn, T., Governato, F., & Stadel, J. 2002, 

MNRAS, 336, 119 
Merritt, D., & Quinlan, G. D. 1998, ApJ, 498, 625 
Moore, B., Kazantzidis, S., Diemand, J., & Stadel, J. 2004, MN- 
RAS, 354, 522 
Naab, T., & Burkert, A. 2003, ApJ, 597, 893 
Naab, T., Jesseit, R., & Burkcrt, A. 2006, MNRAS, 372, 839 
Navarro, J. F., Frcnk, C. S., & White, S. D. M. 1997, ApJ, 490, 
493 

Norman, C. A., May, A., & van Albada, T. S. 1985, ApJ, 296, 20 
Press W. H, Teukolsky S. A., Vetterling W. T. & Flanncry B. P., 

Numerical Recipes in C, Cambridge University Press, New 

York (1993). 
Sideris, I. V. 2006 Phys Rev. E, 73: 066217. 
Stadel, J. G. 2001, Ph.D. Thesis, Univ. of Washington 
Valluri, M., & Merritt, D. 1998, ApJ, 506, 686 
Valluri, M., Vass, I. M., Kazantzidis, S., Kravtsov, A. V., & Bohn, 

C. L. 2007, ApJ in press, astro-ph/0609612 

Warren, M. S., Quinn, P. J., Salmon, J. K., & Zurek, W. H. 1992, 
ApJ, 399, 405 



ACKNOWLEDGMENTS 

All the numerical simulations 

were performed on the zBox2 supercomputer (http://www- 
theorie.physik.unizh.ch/~dpotter/zbox2/) at the University 



6 A. V. Maccid et al. 




5 10 15 5 10 15 5 10 15 5 10 15 

time time time time 



Figure 1. Top row panels: evolution of the axes ratio b/a of two systems, one without central baryonic mass (dotted line) and one 
with central baryonic mass (solid line) for different distances r from the center, (a) r = lA.9kpc, (b) r = 163.4fcpc, (c) r = 311.9fcpe, (d) 
r = AGOAkpc. Bottom row panels: same as top row panels but for axes ratio c/a. (e) r = 14.9fcpc, (f) r = 163.4fcpc, (g) r = 311.9fcpc, (h) 
r = 460.4/cpc. 
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Figure 2. Density contours for the central parts of a dissipationless (top) and dissipational (bottom) halo at the end of the evolution 
(14Gy). There is an obvious rounding of the bottom configuration. 




Figure 3. Complexities (number of important frequencies of motion) for 300 orbits versus mean radial distance from the center. Red 
asterisks represent the dissipationlcss points while black diamonds the dissipational ones. 
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Figure 4. (a) and (b) Poincare sections of the dissipational points recorded at v x = 0. (c) and (d) same for the simulation with a growing 
central baryonic mass. 
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Figure 5. First row panels: x-y, x-z, and y-z projections of an orbit evolving in a halo without a central baryonic mass. Second row 
panels: x-y, x-z, and y-z projections of the same orbit evolving in a halo with a central baryonic mass. Third row panels: Fourier spectra 
of the first (without central mass) orbit. Fourth row panels: Fourier spectra of the second (with central mass) orbit. The Fourier spectra 
in the three dimensions arc broader in the second case. 




Figure 6. Distribution of angles 9 = tan 1 (x/y) explored for the dissipationlcss orbit (black) and the dissipational one (red). 
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Figure 7. Top panels: evolution of the axes ratio b/a of two dissipational systems, one with e = 2.4 (dotted line) and one with e = 0.6 
(solid line) for different distances from the center, (a) r = 14.9fcpc, (b) r = 163.4fcpc, (c) r = 311.9fcpc, (d) r = 460.4fepc. Bottom panels: 
same as top but for axes ratio c/a. (c) r = 14.9fcpc, (f) r = 163.4/cpc, (g) r = 311.9fepc, (h) r = 460.4fcpc. 




Figure 8. Top panels: evolution of the axes ratio b/a of two dissipational systems, one with 500 central masses (dotted line) and one with 
one central mass (solid line) for different distances from the center, (a) r = 14.9fcpc, (b) r = 163.4/cpc, (c) r = 311.9fcpc, (d) r = 460.4fcpc. 
Bottom panels: same as top but for axes ratio c/a. (e) r = 14.9fcpc, (f) r = 163.4fepc, (g) r = 311.9fcpc, (h) r = 460.4fepc. 




Figure 9. Boxcar averaged complexities (number of important frequencies of motion) computed for 300 orbits versus mean radial distance 
from the center. Dark line: dissipational system with only one baryonic mass at the center (e = 0.6). Red line: dissipational system with 
500 baryonic masses at the center (e = 0.6). Blue line: dissipationlcss system (e = 0.6). Cyan line: dissipational system with e = 2.4. 
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Figure 10. Projections of the same orbit for 4 cases: (a)-(c)dissipationlcss system (e = 0.6). (d)-(f)dissipational system with one central 
baryonic mass (e = 0.6). (g)-(i)dissipational system with one central baryonic mass and softening e = 2.4, (j)-(l)dissipational system with 
500 central baryonic masses (e = 0.6). 



